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Abstract 


The thermo-remanent magnetism and the coercive force of the ilmenite- 
hematite solid solution (1-x) Fe,O,;-xFeTiO; have been examined for the range 
0<x<0.7: in the range x>0.7, the solid solution is only paramagnetic at room 
temperature. The specimens are all synthesized ones. Measured quantities are 
J(T), JT), H(T) and H,.T) for the maximum field of 800 Oe. and the TRM 
fit H=2.0 oe <2). The summarized results are as follows: In the parasitically 
ferromagnetic region (x<0.5), /(T) and J,(T) show an increase below the Curie 
point, while H,(7), and H,,(7) decrease monotonously with temperature. Thermo- 
remanent magnetism is strong in comparison with /(7) in this region, the Q-ratio 
amounting to several hundred. In the ferrimagnetic region (0.7 >x>0.5), J/(T) and 
J,(T) show the typical straight decrease with temperature rise. AH, and the Q-ratio 
are quite smaller than for the parasitically ferromagnetic region. Kemarkable re- 
verse TRM like the Haruna-type one appears for the specimens of the range border- 
ing the ferrimagnetic and parasitically ferromagnetic regions, namely the range 
0.6>x>0.45. This indicates that the phenomenon is inherent to this substance and 
is closely related to the fundamental magnetic properties of this series. Further- 
more, it was found that a definite tendency for reverse TRM appears for x++0.1 that 
may be relevant to the natural reverse remanence of Adirondack rocks reported 
by Balsley and Buddington. Stability against the AC demagnetization and the 
field dependence of the thermo-remanent magnetism of this series have also been 


examined systematically. 
1. Introduction 


Ferromagnetism or ferrimagnetism of certain members of the system ilmenite- 
hematite has recently been discovered in the course of the study on the mechanism 
of the reverse thermo-remanent magnetism (reverse TRM) [1]. Since then, a thorough 
examination of magnetic properties of this system has been conducted to establish the 
following facts [2]: when the chemical composition of a member of this system is ex- 
pressed by (1-2) Fe,O,-xFeTiO,, the members in the region 0<*<0.5 are antifer- 
romagnetics with weak parasitic ferromagnetism and those in the 1 region 0.5<%*<0.7 


* Contribution from Division of Geomagnetism and Geoelectricity, Geophysical Institute, Tokyo 
‘University. Series II, No. 69 : 
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are ferrimagnetics at the room temperature. The members with x>0.7 are also ferri- 
magnetics but their Curie point is below the room temperature, whereas the members 
with x+1 are reported as antiferromagnetics with very low Néel point [3]. The com- 
plicated magnetic characteristics of this series are explained in detail in terms of an 
order-disorder transformation of Fe and 7% igns in the crystal structure by Y. Ishikawa 
and S. Akimoto [4]. 

Although the natural occurrence of this series in nature is rarer than that of the 
titano-magnetite series, its importance in rock-magnetism is great because of the various 
peculier characteristics revealed in its thermo-remanent and natural remanent magne- 
tism. For instance, it is precisely members of this series that shows the reverse TRM 
of the Haruna-type [5], [6]. In a previous paper, the author showed that many of 
natural samples belonging to the ferrimagnetic region of this series can acquire the 
capability of producing the reverse TRM by simple heat treatment [7]. Moreover, the 
natural reverse remanent magnetism of the pre-Cambrian rocks of the Adirondack 
Mountains, U.S.A., was suggested to be caused by some yet unknown self-reversal 
mechanism inherent to this series [8]. The pronounced negative geomagnetic anomaly 
at Allard Lake, Ontario, Canada was also suggested to be due to the reverse natural 
magnetization of the vast ilmenite-hematite deposit at that area [9]. The abnormal 
characteristics of this series with respect to remanent magnetism is so marked that it 
has been proposed that any rock containing this series should be excluded from the 
object of the palaeomagnetic studies [10]. 

On the other hand, the natural remanent magnetism of pure kemasite has been 
proved to be useful for palaeomagnetic studies because of its extremely high stability. 
Its thermo-remanent properties has already been studied fairly in detail [11]. 

Considering these facts, it may be noticed that such a systematic examination of 
thermo-remanent properties of this series as has been done for rocks containing 
titano-magnetites by T. Nagata [12] is necessary for further development of rock-magne- 
tism. In the present study, the author has conducted a study to clarify the normal and 

abnormal characteristics of the TRM of this series by using synthesized solid solutions 
- with varying ratio. of ilmenite to hematite. The thermal variation of f the coercive force, 
~ which HS: have important ey to the TRM. 13h has also Dees Pea 
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sand, of which chemical composition had previously been determined [7]. For hematite 
we used two kinds of artificial hematite of different origins: one was synthesized by 
heating of ferric hydroxide which was precipitated by ammonium hydroxide from ferric 
chloride solution, and the other was a commercial guaranteed reagent, both having 
been ascertained to be chemically and crystallographically pure hematite. 

In the following, the series made from the former hematite will be called the 
normal series, and the one from the latter the abnormal series by the reason 
described below. 

The lattice parameters of all the specimens were obtained by a “Norelco” diffracto- 
meter as listed in Tables I and II. The spectrograms of all these specimens show 
exclusively the spectra of rhombohedral crystal structure. Within the limit of the resolv- 
ing power of the “ Norelco,” all specimens were found to be crystallographically homo- 
geneous. Trace of neither the cubic titanomagnetite series nor the orthorhombic pseudo- 


brookite series could be found [2]. 
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Fig. 1 Chemical composition of synthesized Fig. 2 Rhombohedral lattice parameter a,., and 


ilmenite-hematite series, represented on a FeO- Curie point 7, of syntheisized ilmenite-hematite 
TiO,-Fe,O, ternary diagram. series as dependent on chemical composition. 
...normal series, o ....cbnormal series. + ....mormal series, o ....abnormal series. 


The chemical composition obtained by chemical analysis after synthesis of the 
specimens is also shown in Tables I and II and in Fig. 1. Fig. 2. shows the relations 
between the lattice parameter a,, and the Curie point T. obtained by the thermo- 
magnetic analysis and the chemical composition of the specimens. In Figs. 1 and 2, 
the full circles and the hollow ones are for the normal and abnormal series respectively. 
In Fig. 2, it is rather remarkable that, although the lattice parameter a,., and the Curie 
point 7, of the specimens of the normal series are in good agreement with what has 
been determined for the ilmenite-hematite series by T. Nagata and S. Akimoto [2], Grn 
and T,, of the abnormal series considerably deviate from them. Since the values of a,., 
and 7, of the normal series agree with the generally accepted ones, we regard the 
specimens of the xormal series as having the normal characteristics of ilmenite-hematite 
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solid solution and those of the abnormal series as somewhat abnormal. Though the 
observed difference in a@,,, of the two series must have been caused by some subtle 
difference in the nature of the two kinds of hematite mentioned above, it has been not 
possible to clarify the reason any more. The grain size of the specimens used in the 
present study was always kept at about 10/ in mean diameter, so that the differences 


in magnetic properties due to the grain size effect may be regarded as eliminated. 


3. The Method of Measurement 


The magnetic measurements in the present study have been conducted by~a 
ballistic method. The apparatus used is similar to the one reported by us previously 
[14], the major difference from the latter being that the present one is equipped with 
a solenoid coil that can produce a magnetic field as intense as 800 Oe. 

In the first place, the measurements of the thermal variations of magnetization 
J(T), of remanent magnetization /,(7), of coercive force H,(T), and of remanent coercive 
force H,,(T), during heating have been made on each specimen. These measurement 
were carried out by drawing hysteresis loops at successively higher temperatures. 
Definitions of these quantities are as usual and shown in Fig. 3. The remanent coercive 
force, H.,(T), is a coercive force of remanent 
magnetization measured at null magnetic field. 
In the cooling process, the mode of development 
of the total thermo-remanent magnetism has been 
measured in a magnetic field of 2.0 Oe. on each 
specimen. 

For the measurement of the thermal varia- 
tion of coercive force of rocks and ferromagnetic 
minerals; there have recently been reported two 


methods, one by a cathode-ray oscilloscope using 
Fig. 3 Schematic representation of the alternating magnetic field [15], and the other by a 

magnetization curve. magnetic balance [16]. In order to measure other 
quantities simultaneously as in the present study, such a classical ballistic method as 
employed here is most useful. 

The limitation in the present apparatus is that the intensity of the magnetic 
field available by the solenoid coil can not exceed 800 Oe., so that the measured values 
of J(T), JT), HT) and H.,(T) differ from those measured using a magnetic field 
sufficiently intense to saturate the specimen. The results of these experiments will be 
described in §§ 4-7. , ze 

After the series of measurements mentioned above, the stability of the thermo- 
remanent magnetism was examined for each specimen, by means of demagnetization 


by alternating magnetic field. The results are reported in § 8. 


Finally, the dependence of the thermo-remanent magnetism on the intensity of 
the magnetic field applied during cooling was studied for each specimen (§ 9). For 
some specimens, the characteristics of the partial thermo-remanent magnetism which 
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is produed by applying a magnetic field only within some restricted temperature range 
during cooling and of the addition law of thermo-remanent magnetism [12] were also 
examined (§ 7). 

In magnetic measurements of the ilmenite-hematite series, the existence of magne- 
tite impurity must be carefully avoided. Since the magnetism of the titano-magnetite 
series is much stronger than that of the parasitic ferromagnetism of the ilmenite- : 
hematite series, it sometimes is necessary for accurate measurement to reduce the 
amount of the impurity less than 1/100,000 in weight, which is certainly undetected by 
X-ray techniques. Selective dissolution of the titano-magnetite series by hydrochloric 


acid has been utilized for the elimination of it. 


4. Intensities of Magnetization and Remanent Magnetization at Room 
Temperature. [A(7,) and J,(7)).] 


The intensity of magnetization at room temperature /(7)) changes with the 
ilmenite-content (x) but not linearly. The value of /,(T>) varies in the same way as 
J(T,), as observed in Fig. 4, in which circles and triangles represent the values of J(T) 
and J,(T>) for a magnetic field of 800 Oe. The range between ilmenite-hematite is, as 
already mentioned, divided into two parts, namely the parasitically ferromagnetic region 
and the ferrimagnetic region, the border of the two being at about x=0.5. In Fig. 4 
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5. Coercive Force and Remanent Coercive Force at Room Temperature. 


[H-( T») and F.( T»)I. 


In general, the coercive force of the ilmenite-hematite series is remarkably 
stronger than that of the titanogznagnetite series [2], [15]. The variation with chemical 
composition of the coercive force and*the remanent coercive force at room temperature 
is shown in Fig. 5. As observed in this figure, the coercive force decreases with the 
ilmenite content x fairly smoothly while the remanent coercive force shows a remarkable 
increase at x+=0.6. 

As stated already, the maximum magnetic field applied for the present measure- 
ments is only 800 Oe. which is obviously too weak to 
attain the saturation magnetization of the specimens, 
especially in the case of the parasitically ferromagnetic 
specimens. As a result, the values of H, and H., in 
Fig. 5 and Table I are to be regarded as smaller 
than the ultimate coercivities. For example, the 
dependence of the coercive force on the intensity of 


400 800 Hoe the applied field of the specimen x=0.01 (hematite) 
Fig. 6 Example of dependence of js shown in Fig. 6. Extraporating the curve to the 
coercive force H, and remanent ; : ; . 
coercive force H,, on magnetic field Yegion of sufficiently strong magnetic field, the high 
applied. coercivity of this substance reported by, for instance, 


Specimen: hematite (x=0.01) ; 
er R. Chevallier [17] may be expected. 


6. Thermo-Remanent Magnetism. [J a 2.0 Oe. (T,)| 


The intensity of the total thermo-remanent magnetism (TRM) produced in a 
magnetic field of 2.0 Oe. is plotted in Fig. 7 against the chemical composition. In this 
figure, the values of TRM of the specimens of the abnormal series are also plotted by 
hollow cirlces. The most remarkable fact seen in this figure is that the specimens 
x=0.56 and x=0.60 show the “reverse thermo-remanent magnetism” (RTRM). Since 
the chemical composition of these specimens is close to that of the Haruna ferromagne- 
tic ilmenite [7], the reverse TRM of the present artificial specimens may be presumed 
to have the same mechanism as that of the Haruna sample. It was, moreover, shown 
that the synthesized specimens of the ilmenite-hematite +=0.45, 0.48 and 0.51, can 
acquire the same type of reverse TRM by certain heat-treatment. The cross marks 
in Fig. 7 show the reverse TRM of the heat treated specimens. As for the reverse 
TRM, we will discuss in more detail in § 11. Se 

Another remarkable feature of the TRM of the ilmenite-hematite series is its 
great intensity in the parasitically ferromagnetic region. For instance, the intensity 
of the TRM produced in the magnetic field of only 2.0 Oe. are comparable to, and 
gonié times greater than, the intensity of the induced magnetization in the magnetic 
field of 800 Oe. as seen in Tables I and IJ and in Figs. 10~17 and 21~23. This seems 
to be an important characteristic of the ilmenite-hematite series. To show this more 
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Fig. 8 Value of Q-ratio for H=2.0 
' Oe. of ilmenite-hematite series as 
dependent on chemical composition. 
+ ....normal series, 
o ....abnormal series. 


clearly, the Q-ratio for H=2.0 Oe, defined as: 
; TRM 


Qs induced magnetization : 
Fig. 7 Intensity of thermo-remanent is plotted in Fig. 8 against the chemical 
magnetism at room temperature of ae ¥ 
: ilmenite-hematite series as dependent ~ composition. As seen in the figure and 
on chemical composition: TRMs plotted : : RUA gy PA. a 
are the total TRM produced by field- qenles } ang Hy thes Vganags : Bre ree 
cooling in H=2.0 Oe. in the parasitically ferromagnetic region, 
ib eanertg) ceases <0 aaa ~ amounting nearly to 800 for the specimen 
“* x....after heat treatment see text). x=0.17. For titano-magnetite series it has _ 


been known to be quite rare that Q-ratio ; 

exceeds 100 [12]. For the specimens having the Lisp aru of reverse TRM, the 
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ilmenites [7]. In the present cases, the modes of J,(7) and the temperature dependence 
of ultimate coercivities are not known for the magnetic field available is too weak. 
But, at least in the temperature-range near the Curie point, the modes of J wandvot 
H(T) or H.,(T) both seem to be linear as will be seen in Figs. 10~17. 

In Figs. 8 and 9, in addition to the specimens having the Haruna-type reverse 
TRM, specimens x0.1 can be noticed as not being included in the general tendency. 
In the ilmenite-hematite series, the specimen of this chemical composition seems to 
have some particular nature with regard to the thermomagnetic properties, which may, 
in the author’s speculation, be related to the occurrence of the natural reverse remanent 
magnetism of the Adirondack ilmeno-hematites [8] and of the Allard Lake ilmenite- 
hematite deposits [9]. We will treat this in § 12. 


7. Thermal Variation of Magnetization (7), Remanent Magnetization J,(7)), 
Coercive Force H.(T), and Remanent Coercive Force H.,(7), and the Mode 


of Development of Thermo-Remanent Magnetism J 7 (T) 
T., 2.0 Oe 


In the preceding sections, the values measured at room temperature have mainly 
been considered. In this section, the modes of thermal variations of these quantities 
will be examined. 

In Figs. 10~17 the results of the thermo-magnetic measurements on the speci- 
mens of the normal series are shown. In these figures, the symbols are as follows: 

full circles connected by full curve....the magnetization /(7) in the field of 800 Oe. 

full circles connected by dotted curve‘...the remanent magnetization J,(T) pro- 
duced by the field of 800 Oe. 

hollow circles connected by full curve....the coercive force H,(T) after magnetized 

by the field of 800 Oe. 

hollow circles connected by dotted curve....the remanent coercive force H.,(T) 

defined in Fig. 3. 

cross marks connected by broken curve....the mode of the development of the 

total thermo-remanent magnetism Be 20 O ot). 

The general tendency in the thermal variations in /(7), and /,(T) of the para- 
sitically ferromagnetic region of the ilmenite-hematite solid solution series can be 
characterized by their increase with the rise in temperature. This increase, being more 
remarkable for the hematite rich specimens, is considered to be a Hopkinson effect 
which is generally observed in the thermal variation of magnetization of ferromagnetic 
materials when measured in weak magnetic fields. For the parasitically ferromagnetic 
members of the ilmenite-hematite solid solution, having high coercivity, 800 Oe. is 
certainly a weak magnetic field. In Figs. 12~17, the increase in /(T) and J,(T) is seen 

7 to be precisely accompanied by an also remarkable decrease in coercive forces. In this 
series, the coercive forces are seen to decrease with the rise in temperature in a rather 
— manner. In the ferrimagnetic region (Figs. 10 and 11), where the coercive force 
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is comparatively weak, the increase in /(T) with the rise in temperature was not ob- 
served ; J/(7) varied monotonously with the temperature like the natural ferromagnetic 
ilmenites [1]. 

The mechanism for the production of TRM has long been suggested by T. Nagata 
[12] as follows: fine grains of ferromagnetic mineral can acquire intense magnetization 
ina weak magnetic field at some critical temperature where the factors which hinder 
the free movement of magnetic domain walls, such as crystal anisotropy and magneto- 
striction, are diminished, and thus obtained magnetization can be fixed or frozen as the 
specimen is cooled to lower temperature in the*same magnetic field. Néel put such 
idea of freezing of magnetization by weak-magnetic-field-cooling into more explicit forms 
[13]. In Néel’s theories, the blocking of magnetic domains during field cooling from 
above the Curie point is presumed to take place at a temperature where the coercive 
force of the sample becomes equal to certain critical values; in one model the critical 
value is the externally applied magnetic field and in another model it is the thermal 
agitation field. In any case, it will be of significance to determine the relation between 
the critical temperature of thermo-remanent 
magnetism and the coercive force. Fig. 18 
shows, for instance, the partial thermo- 
remanent magnetism characteristics and the 
coercive force of a specimen x=0.01 (hema- 
tite). 

In this figure, it is clearly observed 


200 that the effective temperature range for the 


| 


production of TRM is quite narrow and 
close to the Curie point; over 90% of the 
total TRM is produced in a temparature 


range of 20° including the Curie point. It 
is also found that in this range, the coercive 
force becomes of the comparable magnitude 
with the external magnetic field that is 10 
Oe. in the present case. The addition law 
of TRM, [12], was also ascertained to be 


100 


Total TRM 


ete Cee 700°C valid for this specimen as the figure shows. 
Fig. 18 Example of the spectrum of partial 
; T-20 
thermo-remanent magnetism J7- “49 (¢ (To) 8. Stability of Thermo-Remanent 
and thermal variation of coercive forces near Magnetism of Tim enite-Hematite Series 


the Curie point: specimen. .hematite (x=0.01). 


Vertical lines at the left side of the figure against the Demagnetization by 


represent the summation of the partial TRMs Alternating Magnetic Field. 
(full line), and the total TRM (broken line) 
respectively, indicating the validity of the addv- Fig. 19 shows the demagnetization 


ion I TRM. he i i i 
Re hoeeah curves of TRM of the ilmenite-hematite 


series by alternating magnetic field. The TRM tested has been produced by a total 
magnetic field cooling in a magnetic field of 2.0 Oe. The ordinate, therefore, represents 
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Jr) #/Ure)#-9 in the conventional notation 
[10]. In this figure, it is observed that the 
normal TRM of the ferrimagnetic samples 
(x=0.71 for instance) is not very stable 
against AC demagnetization, whereas the 
normal TRM of the parasitically ferromagne- 
tic samples (x<0.3) and the reverse TRM are 
nearly completely stable. 

In comparison with the similar curves 
for rocks containing titanomagnetite series . 
reported by us previously [10], the stability 
of TRM of the parasitically ferromagnetic 
Fig. 19 Stability of thermo-remanent magne- samples revealed in Fig. 19 is very high. 
tism a 2.0 Oe (Lo) of ilmenite-hematite This high stability of TRM is a natural con- 
series against demagnetization by alternating sequence of the high coercive force of these 
magnetic field H. samples. 

It is a very remarkable fact that the reverse TRM of the ferrimagnetic specimen 
is as stable as the TRM of the parasitically ferromagnetic ones. In Fig. 19 the double 
circles show the reverse TRM of the specimen x=0.48 after heat-treatment (explanation 
in § 12). 


9. Field Dependence of TRM of Ilmenite-Hematite Series. 


Saturation phenomenon of TRM with the increase in the externally applied mupie. 
tic field has been recognized for rocks and _ . 
bricks by various authors. T. Nagata empiri: —§  emuX¢ —~ omy %¢ > 
cally formulated this phenomenon by, [12], . 


7 eates inh CH. } (a). 


where kisa constant, and L. Néel provided 
ary theoretical basis for it by a fine particle 
0 lel 12). Recently i Néel | proposed alter- 

expression for this phenomenon by 
eo mite sf vil 
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90 Oe., the apparently perfect saturation of TRM for the field greater than some 100 
Oe. shown in the figure is quite remarkable. 

In this figure, it is also to be noticed that the specimens x=0.60 and x2+0.1 are 
again of somewhat different nature from other specimens: their saturation effect is 
not as perfect as that of other specimens. The double circles again stand for the 
reverse TRM of the specimen x=0.48 after it acquired the reverse TRM characteristics 
by heat-treatment. It will especially be noteworthy that the specimen x=0.48 can 
retain its reverse TRM characteristics even in the field as intense as 500 Oe., because 
this fact means that the effective magnetic field orienting the TRM of this sample in 
the reverse direction must overcome the normally directed external magnetic field of this 
strength. The bearings of this fact will be treated in § 11. 


10. Summarized Results of the Specimens of the Abnormal Series. 


So far, we have mainly referred to the experimental results obtained for the 
specimens of the zormal series. In this section, the corresponding results on specimens 


of the abnormal series defined in § 2 will be reported briefly. Figs. 21~23, are the 
examples of the results of the thermo-magnetic measurements, where the symbols are 
identical with those in Figs. 10~17. Although the increase in J/(T) and in J,(T) with 
rise in temperature is similar to the case of the xormal series, the increase in the 
coercive forces of the parasitically ferromagnetic specimens illustrated in Fig. 22 is in 
complete contrast to the case of the zormal series. Among the specimens with «<0.5 
of the abnomal series, only specimen x=0.13 (Fig. 23) did not show such increase in 
HAT) and H,,(T). We have not yet-found any plausible explanation for this increase 
of H, and H..,. The values of J/(T), J(To), H,(To) H.»(To) and the @-ratio are listed 
in Table II. Probably due to the above queer nature of the coercive force, such a simple 
relation between TRM and coercive forces as in Fig. 9 does not exist for this series. Fig. 
24 shows the field dependence of TRM of the abnormal series. Nearly perfect saturation 
phenomenon of TRM is observed also for this series. The slower saturation for the 
specimens x=0.56 (RTRM) and x+<0.1 is similar to the corresponding cases in Fig. 20. 
As was mentioned at the end of § 6, the specimens x+=0.6 (reverse TRM) and x+0.1 
showed the abnormal behaviours in TRM in the normal series. It is now clear that 
the same is true also for the abnormal series. (see also Fig. 7) 


11. Reverse TRM Characteristics of Ilmenite-Hematite Series. 


As stated already, some members of the ilmenite-hematite series show the reverse 
thermo-remanent magnetism. (see Table I.) In Fig. 25, the directions of TRM of the 
present specimens are shown by N (normal) and R (reverse) attached to the points 
representing the chemical composition of each specimens. As for the symbol R’, ex- 
planation will be given in § 12. The specimens with x=0.60 and 0.56 show a 
marked reverse TRM of which mode of development, shown in Figs. 11 and- 
21, is just similar to that of the reverse TRM of the natural ferromagnetic ilmenite 
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Fig. 21 Fig. 22 Fig. 23 


Figs. 21-23 Examples of thermal variation of magnetization /, remanent magnetizaton /,., 
coercive force H,, remanent coercive force H,, of the abnormal ilmenite-hematite series for 
H=800 Oe. and their mode of development of thermo-remanent magnetism produced in 
-H=2.00e. Symbols are same as in Figs. 10-17. 
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Sample in the Haruna pumice [5]. Since the chemical composition of the natural 
ferromagnetic ilmenites which have the Haruna-type reverse TRM are known to be 
x=0.6~0.5, (double circle in Fig.,25 represents, for instance, the Haruna sample), it is 
quite probable that the reverse TRM of the specimens synthesized in the present study 
and the reverse TRM of the Hardna-type have a common origin. Since the specimens 
examined in the present study have been, during the process of synthesis, quenched 
from 1200°C, which definitely is above the top of the solvus curve of the ilmenite- 
hematite series [18], the two-phase hypothesis postulated hitherto [5] to explain the 
various characteristics of the Haruna-type reverse TRM must be reconsidered. 

Moreover, the field dependence of the reverse TRM of the present specimens 
revealed in Fig. 20 indicates that the reverse TRM is still acquired by cooling in 
magnetic field as intense as 500 Oe. It can be readily shown that such an intense 
reverse magnetic field can hardly be accounted for by demagnetizing field in the 
conventional two-phase hypothesis [5].* 

It seems likely that such intense reverse effective magnetic field may be caused 
by some coupling mechanism related to the Weiss-Heisenberg exchange forces. 

In Fig. 7 or in Fig. 25, it is remarkable that the specimens showing the reverse 
TRM are concentrated in a limited region near the boundary between the ferrima- 
gnetic and parasitically ferromagnetic regions of the ilmenite-hematite series: namely 
0.45<x*<0.6. The magnetism of the specimens of this region has been noticed to 
be very sensitive to heat treatment [19], and on the basis of this sensitiveness on 
heat treatment, the production of ferromagnetism of the ilmenite-hematite series has 
been attributed to an ordering of 77 and Fe ions in the crystal lattice [4]. As mentioned 
in § 6, the characteristics of the reverse TRM has also been found to be sensitive to 
heat treatment. The best example of this is as follows. 

The specimen denoted by a cross mark in Fig. 25 is a synthesized one from 
TiO., Fe,O; and Fe in the same method as others and its chemical composition and the 


lattice parameters are: 


FeO: 31.20 % in mol 0, =3ATIAA 
Fe,0; 35.19 ” Arn 55°05" 
AiOssn 33:61 is a t= 048 


This specimen (x=0.48) was parasitically ferromagnetic and had only the normal TRM 
characteristics immediately after the synthesis. But, after being kept at 1000°C for 8 
hours in an evacuated sealed quartz tube and slowly cooled, (this heat treatment is known 
to be effective in promoting the ordering for production of the ferrimagnetism of the 
specimen of this composition [19]), it became ferrimagnetic and acquired the reverse 
TRM characteristics. The reverse TRM of the specimens with x=0.48 in the normal 
series and x=0.51 and 0.45 in the abnormal series have also been acquired by the same 


* For instance, the possible maximum demagnetizing field expected from an alternating parallel 
structure is Hg=2zJ, where J is the saturation magnetization of the component with higher 
‘Curie point. For this value to exceed 500 Oe., J must exceed 80 emu/cc whereas the saturation 
magnetization of our present samples is much smaller. 
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heat treatment, their original TRM being normal as shown in Figs. 7 and 12 and 
Tables I and II. It may be, therefore, well concluded that the production of the reverse 
TRM is closely related to the ionic ordering essential to the production of the ferri- 
magnetism of the ilmenite-hematite series. Since it is probable that the specimens 
in the boundary region between the ferrimagnetic and parasitically ferromagnetic 
regions contain both the ordered and disordered phases, the assumed Weiss-Heisenberg 
exchange interaction may be coupling these two phases, as in the case of Co and CoO 
reported by W.H. Meiklejohn and C.P. Bean [20]. 

Emphasis should, now, be placed on the fact that even purely artificial specimens 
show -the reverse TRM. We shall discuss more fully on the mechanism of the reverse 


TRM in a forthcoming paper separately. 


12. Imperfect Reverse TRM Characteristics of Ilmenite-Hematite Series. 


The specimens with x+0.1 in both the normal and abormal series show an 
anomalous characteristic with respect to the mode of the development of TRM (cross 
marks in Figs. 16 and 23.). In an ordinary mode of the development of TRM, the 
TRM increases monotonoulsy with the decrease in temperature. But, the curves in Figs. 
16 and 23 show that the TRM decreases from about 200°C. Although the TRM of 
these specimens at the room temperature is still normally directed, it definitely has a 
tendency for reverse TRM: such TRM will be called “imperfect reverse TRM ” here- 
after. The imperfect reverse TRM is denoted by R’ in Fig. 25. 

That this tendency for reverse TRM is 
not caused by the presence of the component 


x= 0.13 


having the Haruna-type reverse TRM as 
impurity will be proved by the following ex- 
periment: in Fig. 26, the full and the hollow bits 100 ce 
curves represent the pee of the Seley 
of the partial TRMs Feri 100 Oe and Tere 

100 Oe respectively of the specimen x«=0.13 

(abnormal series). From these curves, it can 
be stated that the decreaseof TRM, chara- 


cteristic to the imperfect reverse TRM, takes 


ae 200 400 600°C 
place by the partial field-cooling in the high pee es sf ——; 


temperature range such as between 570°C. Fig. 26 Modes of development of partial 
500°C, and not in the range between 400°C. *hetmo-remanent magnetisms Jer0°, 100 Oe. 
250°C. The lower temperature range is nd Licah. 100 Oe, Of the specimen x= =0. 13. 
known to be the effective range for produc- 
tion of the-Haruna-type reverse TRM [5]. 

From this experiment, it is clear that for the production of the imperfect séverte 
TRM, the normally directed TRM of the phase having a Curie point at about 560°C 


is nesessary, and that the phase with lower Curie point is magnetized reversely by some 


Thermo-Remanent Magnetism and Coercive Force of the Ilmenite-Hematite Series. 77 


interaction with the already produced normal TRM of the high Curie point phase, 
which is the major constituent of the specimen x+0.1. 

An alternative explanation for the phenomenon may be possible if we assume 
that the thermal variation of the saturation magnetization of the ilmenite-hematite 
series is Néel’s P-type [21] at thé composition x-+0.1. 

According to J.R. Balsley and A.F. Buddington [8], the magnetic minerals in the 
reversely magnetized rocks at Adirondacks are hematite rich ilmenite-hematite solid 
solutions with exsolution lamellae of ilmenite. The present author, therefore, considers 
that the natural reverse remanent magnetization of the Adirondack rocks may be 
caused by the similar mechanism as that of the imperfect reverse TRM of the present 
case, in the Adirondack rocks the extremely slow cooling in the process of metamor- 
phism possibly having made the mechanism work so effectively as to bring about the 
perfect reversal of the remanent magnetism. 


13. Conclusion 


The thermo-remanent magnetism and the thermal variation of coercive force of 
the ilmenite-hematite solution series (1—x) Fe,O;-xTiFeO; have been investigated on 
synthesized specimens for the range 0<x*<0.7; the range x>0.7 has not been examined 
because the solid solution is known to be only paramagnetic in this range at room 
temperature. The range 0<x*<0.7 is known to be divided into two ranges: the ferri- 
magnetic range 0.7>x>0.5 and the parasitically ferromagnetic range 0.5>«>0 accord- 
ing to the intensity of magnetization. The magnetic characteristics investigated in the 
present study, such as the intensity and the thermal variation’of the induced and 
remanent magnetization, of coercive forces, and of the thermo-remanent magnetism, 
were also found to be quite different in both ranges. In ccmparison with the ferri- 
magnetic range, the parasitically ferromagnetic range is characterized by the weak 
magnetization, high Curie point, high coercive force, high @Q-ratio, and high stability 
of TRM against AC demagnetization. 

With respect to the characteristics of the thermo-remanent magnetism, the border 
range (0.45<x<0.6) between the ferrimagnetic and parasitically ferromagnetic ranges 
is quite conspicuous for the occurrence of the reverse thermo-remanent magnetism. 

The specimens with x2+0.1 are also remarkable because it shows a definite tendency 
for reverse TRM. 

The ilmenite-hematite solid solution is, thus, quite important in rock-magnetism 
because of these various abnormal thermo-remanent characteristics of which origin is, 
from various experimental results, considered to be closely related with the very funda- 
mental nature of the magnetism of this series itself. Further studies are under progress 


- 
a 


along this line. 
In concluding, the author would like to express his deepest gratitude to Prof. 


T. Nagata for his kind directions throughout the course of this study. He also wishes 
to thank Dr. T. Rikitake for his encouragement, to Dr. S. Akimoto, Dr. S. lida and 
Mr. Y. Ishikawa of Tokyo University for their valuable discussions and to Mr. M. 


78 S. Uvyrepa 


Ozima, now at University of Toronto, for the construction of the apparatus for the | 
measurement. For the chemical analysis of the synthesized specimens, the writer is 
wholly indebted to Dr. K. Katsura of Tokyo Institute of Technology. 
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Abstract 


Observations on the correlation of solar activity and cosmic radiation are 
reported. The relation of these results to models involving solar corpuscular 


streams is briefly conidered. 


Introduction 


Investigations of the time variations of cosmic radiation have established a 
connection between some aspects of solar activity and cosmic radiation; in particular, 
it has been found that a few large solar flares are accompanied by large but brief 
intensity increases on the earth and on some occasions are followed by decreases, 
together with magnetic storms.! In addition to these variations, evidence has been 
obtained which indicates that there are variations associated with the 27-day period 
of solar rotation? as well as the 1l-year period of the solar cycle.® 

Since the solar flare component as well as the main component, presumably of 
galactic origin, approach the earth through the solar system, both components of 
cosmic radiation are subject to the influence of interplanetary magnetic fields. Thus, 
the electromagnetic conditions in the solar system are expected to be reflected in such 
experimental observations as the rate of rise and fall of flare radiation, the 27-day and 
1l-year variations. 

In this connection, recent investigations‘-§ have examined questions relating to 
the source, nature and role of interplanetary magnetic fields in cosmic ray effects. 
The fields considered in these discussions have their origin in the galaxy as a whole 
or on the sun and are turbulent in nature, due to either the pressure of solar 
corpuscular radiation or the hydromagnetic transport of solar fields. Field configura- 
tions which have been treated include 1) a field-free cavity,4 swept out by solar 
streams, with walls made up of lines of force of the galactic field, 2) turbulent solar 
streams of magnetized gas clouds*-? and 3) spherical shells®* of tangled or disordered 
fields centered on the sun. Whichever model is finally adopted, it must satisfy certain 
specific cosmic ray requirements: 1) allow relatively direct orbits connecting the sun 
and the earth so as to provide a means whereby the rise of flare radiation at the 
earth may be fairly rapid, yet 2) provide some. means to explain the slow decay of 


* Now. at the University of California, Berkeley, California. 
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flare radiation and further 3) provide modulating mechanisms to account for the 27-day 
and ll-yeay variations. Other requirements involve astrophysical considerations: the 
model must avoid 1) degrees of symmetry which cannot be maintained in nature and 
2) excessive drains on solar resources in the form of corpuscular emission and 
magnetic fields. 

While the details remain to be worked out, the model®-? involving magnetized 
streams emanating from the equatorial regions of the sun appears to offer promise. 
For this configuration the drain on solar resources is reduced considerably and the 
geometry is consistent with the known features of solar activity. In addition, some 
progress has been made with this model toward understanding flare effects by consider- 
ing the diffusion of particles through magnetized gas clouds and modulation effects by 
considering the capture and support of magnetized gas clouds by the earth’s gravita- 
tional and magnetic fields. 

It is the purpose of this paper to present experimental observations during a 
period of intense solar activity which may have bearing on these questions; in 
particular, to consider the relation of solar activity and slow cosmic ray variations as 


shown in the period surrounding the 1956 flare outburst.® 


Observations 
A. Recent Event 


Forbush? has published the results of a long-term study of the variations of 
cosmic ray intensity registered by Compton-Bennett ionization chambers. The data 
obtained during the period 1937-1952 show that there was a variation of cosmic ray 
intensity over the 1l-year solar cycle, annual means of cosmic ray intensity displaying 
a negative correlation with sunspot number. In this study the possible influence of 
short term magnetic disturbances was examined by comparing the variation of annual 
means for all days, international magnetic quiet and disturbed days. The results of 
this comparison indicate that variations of intensity with sunspot number persist for 
periods longer than those characteristic of transient magnetic disturbances. 

Kodama and Murakami’? have published results of neutron intensity observations 
at Nagoya for the interval. September 1954 to August 1955. Their observations, 
although obtained near solar minimum, also indicate a negative correlation of cosmic 
ray intensity with sunspot number. 

Since it has been established that a linear relationship connects long-term 
averages of sunspot numbers and areas, a correlation similar to that with snnspot 
number is expected to exist between cosmic ray intensity and sunspot area. In addi- 
tion, it is to be expected that the size of sunspots would be more closely related to 
solar-terrestrial effects than number, especially during periods of intense solar activity — 
' lasting times short compared to the length of the solar cycle. In view of these 
points, it is of interest to consider the solar cosmic ray outburst of February 23, 1956 
which was preceded by very intense solar activity. 
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Observations of the cosmic ray intensity variations surrounding this event, !? 
obtained with a neutron monitor, indicate a close time-association between the daily 
averages of cosmic ray intensity, corrected for barometric fluctuations, and the disk 
passage of active regions on the sun. In addition, examination of sunspot data!® for 
this period, as shown in Fig. 1/ indicates that the cosmic ray intensity was negatively 
correlated with sunspot area, at least for the months of January and February. 


Fig. 1 Cosmic Ray Intensity 
(nucleonic component) and 
Sunspot Area for the Period 
January-April 1956 (Area ex- 
pressed in millionths of the 
sun’s hemisphere and _ cor- 
rected for foreshortening). 


counts /MiNUTE 


SUNSPOT AREA 


The role of transient magnetic disturbances is more difficult to determine in a 
short-term study such as with this event. One possible approach is to indicate the 
magnetically quiet and disturbed days!‘ during the period and examine their distribu- 
tion relative to the cosmic ray variations. These days are given in Table. 1 and are 
indicated by Q (quiet) and D (disturbed). Inspection of the data in Fig. 1 shows that 
for the months of January and February the intensity variations were relatively slow, 
both during the decreasing and increasing phases. Further, the distribution of Q and 


D-days given in Table. 1 does not show Table 1. Magnetically quiet and disturbed 
any unusual trends which might suggest days.4 

that ue eva were due to \the Month eis eee: ae 
superposition of a number of transient ie 
magnetic storm effects; rather, it is seen s | 1 

that Q-days are found to occur while the January ae 2 
intensity was decreasing and D-days while 26 28 

it was increasing. Z " 

For the months of March and April, February f = 
following the cessation of intense solar 14 29 
activity, the cosmic ray intensity varia- ee z 3 
tions were of a different character; in- March 9 29 
stead of a slow decline and recovery, us - 

lasting about 14 days, the intensity ny eo" 21 
decreased rapidly over two or three days | April 15 27 
and then recovered slowly, over a period i a 


of 45 days, with brief variations super- 
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imposed. These variations are probably related to the next disk passage of the active 
regions responsible for the variations in January and February. Toward the end of 
April an intense geomagnetic storm occured; this storm, in contrast to the magnetic 
disturbances which occured earlier,'? produced the familiar rapid decrease and recovery 
in intensity. : 

Just as the intensity variations for March and April differ from those in January 
and February, so does the correlation between intensity and sunspot area. Thus, 
inspection of Fig. 1 shows the sunspot area considerably above the minima for January 
and February but in poor correlation with the intensity variations. 
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whereas the March-April data show considerable scatter and deviation below the 
correlation line. In these plots the magnetically quiet (Q) and disturbed (D) days are 
represented by circles and squares, respectively; it is seen that their distribution does 
not show any pronounced trends. 

Cosmic ray intensity observations obtained with an ion chamber located at the 
geomagnetic equator!® (Huancayo, Peru) also show the intensity variations described 
above. A presentation of the 1956 ion chamber data, similar to that given for the 
neutron monitor data, is shown in Fig. 4. While the ion chamber observations 
resemble rather closely those obtained with the neutron monitor, it is seen that the 
amplitude of the variations is considerably smaller. This is due to the location of the 
detector at a low geomagnetic latitude and the fact that the ion chamber responds 
mainly to the mu-meson component which has a higher mean primary energy than 
the nucleonic component. 


be | Fig. 4 Cosmic Ray Intensity 

(meson component) and Sun- 

spot Area (expressed in mil- 

| lionths of the sun’s hemi- 

| sphere) for the Months of 
January and February 1956. 
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B. Previous Events 

Four solar cosmic ray outburst of the type observed on February 23, 1956 have 
been observed earlier: February 28, 1942, March 7, 1942, July 26, 1946 and Novemker 
19, 1949. With the exception of the November 19, 1949 event,’® all previous events 
have been observed with Geiger counter telescopes or ion chambers. As mentioned 
above, these devices respond mainly to the mu-meson component for which the average 
primary energy, at a given geomagnetic latitude, is considerably greater than that for 
the nucleonic component. It has been found in the case of these events that flares 
produced increases in intensity of the order of 30% in such detectors; in addition, the 
- increases were followed by a slow decay lasting seveal hours after the flares. However, 
in contrast to the 1956 event, no increases were observed at the geomagnetic equator. 
Intensity decreases due to magnetic disturbances following some of the flares, on the 
other hand, were world-wide, being observed at middle as well as low geomagnetic 
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latitudes. 
It would be of interest to examine the data’ for these earlier events for modula- 


tion effects similar to those observed prior to the 1956 event. Unfortunately, the ion 
chamber data!® and the solar activity data,!’ while not inconsistent with solar-controlled 
downward modulation, do not show clear-cut intensity decreases for these events. 
Considering the amplitude of the variations at Huancayo in the case of the 1956 event 
and the fact that the solar activity preceding the earlier events was less intense and 
of shorter duration (expressing sunspot area in millionths of the visible hemispher, 
Amar~2500 for the 1942 events, Amar~5000 for the 1946 event and Amar~3500 for the 


1949 event), this result is not unexpected. 


Discussion 


One feature common to all the flare events is the slow decay of post-flare radia- 
tion. At the present time there seems to be two possible explanations of this feature: 
either a continuous but gradually decreasing emission of flare particles, extending 
beyond the duration of the visible aspects of the flares, or the trapping and leakage 
of flare particles out of an interplanetary field. Of the five outbursts, the 1956 event 
has been studied in the greatest detail. Experimental data®!*® on the distribution of. 
the post-flare radiation for this event indicates that it had a rather high degree of 
isotropy and thus argues against the continuous emission of particles by the sun. 

Whether the interplanetary magnetic field, required by the above discussion, has 
a relatively permanent existence or is occasionally transported out from the sun by 
streams resulting from enhanced solar activity, is not clear from observations of earlier 
flare events. However, the striking correlation of cosmic ray intensity with sunspot 
area and disk passage of active regions in January and February strongly suggests 
that, in this case, these decreases find their explanation in modulation effects produced 
by magnetized solar streams. These observations seem consistent with a geocentric 
model involving the interaction of the earth’s magnetic field with a broad solar stream, 
perhaps along the lines indicated by Parker® where modulation effects are ascribed to 
the capture and support of additional magnetized clouds sent out from the sun by the 
earth’s gravitational and magnetic fields. Within the framework of this interpretation 
the sunspot area of the visible hemisphere would give a measure of the density of the 
streams emanating from the sun from which clouds are captured by the earth. 
However, the duration of the decreases would be related to the details of the capture 
and retention of clouds from the streams. In view of the fact that the intensity 
variations follow the solar activity rather closely, it would appear that the clouds 
captured by the earth’s gravitational field are not bound very strongly by the earth. 
In developing an explanation of this feature, it would seem that capture of clouds 


would play a dominant role while the solar stream is overtaking and enveloping the 


earth while the loss of clouds through collisions with clouds in the stream would 
become RACE ae ee the earth has passed through the more central portions of os 
stream. . 
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Abstract 


In Parts I and II ionospheric wind systems have been deduced from the 
geomagnetic Sg variations for the mean solstice. In order to compare our results 
with observation, however, it is still necessary to obtain the wind systems for 
the solstitial seasons. In this paper wind systems in summer and winter hemis- 
pheres are deduced from the Sy data. It is shown that the diurna] wind component 
has greater magnitude in summer than in winter; whereas the semi-diurnal 
component has smaller magnitude in summer than in winter over most of the 
hemisphere, and that in summer hemisphere the diurnal wind component is 
greater than the semi-diurnal one; whereas in winter hemisphere the former is 


more or less smaller than the latter. A comparison between our results and 
observed ionospheric winds is made. 


1. Introduction 


In Part I (by this author [1]) and Part II (by S. Kato [2]) of this report we have 
obtained the horizontal wind systems in the ionospheric EF region from the dynamo 
theory of the geomagnetic S, variation, and it has been found that the main component 
of winds required to produce the S, variations is not semi-diurnal, but diurnal. M. 
Hirono and T. Kitamura [3] have also treated this problem on some different assump- 
tions, and obtained the results similar to ours. Since all these works [1,2,3] are based 
on geomagnetic data for the mean solstice, the results seem to represent a mean state 
of summer and winter. In order to compare the calculated results with observation, 
it is, therefore, still necessary to obtain the wind systems for each of both solstitial 
seasons. In the present paper (Part III) the wind systems in Northern summer and 
winter (i.e. Southern winter and summer, respectively) are deduced from the geomagne- 
tic data by a similar method to in Part I, on the assumption, for simplicity, that the 
motion of the air in regions concerned is horizontal and irrotational. Further extension 
of Part III, taking into account the effect of the Coriolis force, will be made in Part 
IV by Kato [4]. . 

The data used are the horizontal intensity of the geomagnetic S, variations at 62 
stations over the world in summer and winter of the Second Polar Year, 1932-33. The 
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harmonic coefficients of the North (4X) and East (4Y) components of S¢ are plotted 


in Fig. 1. The full (for Northern’ summer) and dotted (for Northern winter) lines in 
this figure show a smooth distribution of the coefficients derived from the potential 
of the S, field for the mean solstice [5] and solstitial difference [6] (see Table i), ore 
the assumption that the mean solstitial field of So is Symmetrical, and that the solstitial- 
difference field is anti-symmetrical With respect to the equatorial plane. Since the 
spherical harmonic coefficients shown in Table 1 are obtained from the S, data in the 
interzonal region (between 60°N and 60°S), the high-latitude curves in Fig. 1 are an 
extraporation of S, in middle and low latitudes. 
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Fig. 1 Latitude-distributions of the harmonic coefficients for the geomagnetic North (AX) 

and East (4Y) components of the geomagnetic Sg variations in surares and winter. The 

full (for Northern summer) and dotted (for Northern winter) lines show a smooth 
distribution obtained from the spherical harmonic coefficients listed in Table 1, 
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Fig. 2 Latitude-distributions of the constant term Cp. 


Table 1. Spherical harmonic coefficients for the mean solstitial and solstitial-difference fields 
of Sg in the interzonal region (between 60°N and 60°S) analyzed in the form 


VE > (an cos mt +b” sinmt) P?? (cos 6), 
using the Second Polar Year data. (Unit=7) 
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‘ Done 4 Dero 
Dyy f=4 D5 yvo (1+ 2.00 cos x + 1.46 cos? x), 
Dey ‘ Deyo 


where . 
x= solar zenith angle. 
The total field for the Northern summer (i.e. Southern winter) is shown in Fig. 3. 
Thus the differential equation (1) for the electrostatic potential S can be solved. 
The boundary conditions satisfied by S are: S must vanish at the poles, and the static 
field derived from S is equal to the total field at the equator, because the dynamo field 
vanishes there. . 
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Emu. 
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Fig. 3 Latitude-distributions of the 
harmonic coefficients for the total 
electric field components E,, Ey in 
Northern summer. 
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S (0, 2)=33 (Spa(O)-COS MA +Sy,(8) Sin mA}, 
then (1) is reduced to 


0s 1 ds m* 
mnt : ma — ____§) = fing (9), 
00? sin@ cosd 00 sin’ @ f | 


OS ng 1 OSmp bare 0). 
ag? sind cos) 00. si in?8 Soa Tend 


(3) 


The general solution of (3) consists of a particular solution and complementary 
functions. A particular solution sm» is obtained by means of numerical integration. 
The homogeneous equation of (3) is reduced to the hypergeometric equation by the 
same substitutions as in Part I. One of linearly independent solutions of this equation 
tends to infinity as 0>0, hence only the following solution is adopted: 


ant pasinnor( 2, Sila m+1; sin*@) 
2 2 , 
where F(a, b;c;z) stands for the hypergeometric function. The general solution is, 


therefore, given by 
Sin = Simpy tz G Smileys 


where C is an integration constant. Since both particular and complementary solutions 
satisfy the conditions at the poles, the value of integration constant C can be deter- 
mined by the condition at the equator, consequently the unique solution is obtained. — 
Further detailed descriptions of the method of obtaining the solution will be given in 

: Seats veby Kato [4]. 


4. Results = 
The static field H, is derived from the potential S ee follows: al a 
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is towards the North (South) at 1500 hr in summer in the Northern (Southern) hemis- 
phere, and at 1100 hr in winter in’the Northern (Southern) hemisphere. 3) The magni- 


tude of the rotating vectors is about 40 m/sec in summer, and about 15 m/sec in winter, 
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Fig. 4. Distributions of the wind velocity in Northern summer for each of the diurnal (left) 
and semi-diurnal (right) components. 
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_. Fig. 5 Distributions of the wind velocity in Northern winter for each of the diurnal (left) 
and semi-diurnal (right) components. 
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for both hemispheres. 
(b) Semi-diurnal wind component 

1) The 12-hourly harmonics represent a clockwise rotation in the Northern hemis- 
phere, and an anticlockwise rotation in the Southern hemisphere. 2) The movement 
is towards the North (South) at 1130 and 2330 hr in summer in the Northern (Southern) 
hemisphere, and at 0230 and 1430 hr in winter in the Northern (Southern) hemisphere. 
3) The magnitude of the rotating vectors is about 10 m/sec in summer, and about 
20 m/sec in winter, for both hemispheres. 
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clockwise in the Southern hemisphere. This direction of the rotation and the magnitude 
of the rotating vectors are in general agreement with the observed results [8] (for both 
harmonics) [9, 10] (for semi-diurnal). However, the phase of the rotation does not agree 
with the observed one. The observed phase is such that the diurnal wind component 
blows towards the South near @600 hr in the Southern hemisphere [8], and that the 
semi-diurnal wind component blows towards the North (South) near 0600 and 1800 hr 
in the Northern (Southern) hemisphere [8, 10]. The difference between the calculated 
and observed phases may be accounted for by the fact that the winds at altitudes of 
80-100 km which are probably lower than the dynamo region are measured by the 
radio-echo meteor technique. 

The diurnal wind vectors have greater magnitude in summer than in winter, 
whereas the semi-diurnal wind vectors have smaller magnitude in summer than in 
winter over most of the hemisphere. And, in summer hemisphere the diurnal com- 
ponent of winds is greater than the semi-diurnal one, whereas in winter hemisphere 
the former is more or less smaller than the latter. This situation seems to be consis- 
tent with the results observed by L. Harang and K. Pedersen [11]. For the sake of 
exact comparison, however, further systematic measurements of ionospheric winds, 
especially in the dynamo region (100-130 km), are still necessary. 
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Geomagnetic Storms Ill. Auroral Latitudes 
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Geophysical Institute, Kyoto University 
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Abstract 


A study is made, based on the electron drift theory, of ionospheric F2 dis- 
turbances in the auroal latitudes associated with geomagnetic storms. The process 
of the study is analogous to that in middle and lower latitudes as shown in the 
previous papers J and II. The results show that F2 disturbances in the auroral 
latitudes, as well as those in lower ones, are ascribed to the effect of the vertical 
drift of the electron caused by electric field deduced from the geomagnetic dis- 


turbance-daily variation. 
1. Introduction 


Ionospheric F2 disturbances in the auroral latitudes associated with the geomagnetic 
storm have not yet been fully studied statistically. One of the reasons is that on a 
magnetically disturbed day a full observation on the F2 region is impossible owing to 
absorption of the radio wave or an obscuration due to the sporadic E region, hence 
the disturbance variations of the fF2 and h,F2 (critical frequency and the height of 
the maximum electron density of the F2 region) are difficult to obtain. In spite of 
these circumstances, some results for variation of the AF 2 are obtained by Nagata [1], 
Martyn [2] for an average state and by Meek [3], Obayashi [4] for individual ones. 
A variation of the h’F2 (virtual height) is obtained by Martyn [2], but, the result 
seems to be doubtful, since retardation of the radio wave was not taken into cosidera- 
tion. On the other hand, theoretical studies have been made by Nagata [1] Meek [5] 
(thermal expansion theory) and Martyn [2], [6], Maeda [7] (ionization drift theory). 
But their treatments are confined mainly to supposition or rough estimations of the 
order of variation, and detailed calculations have not yet been made. 

We have shown in two previous papers [8], [9] (hereinafter referred to as I and II) 
with the same title as the present that the F2 disturbances in the latitudes below the 
auroral can be well explained as an effect of the vertical drift of the electron due to 
the electric field deduced from the disturbance-daily (Ds) variation of the geomagnetic 
field. In this paper we further attempt to explain the auroral F2 disturbances as the 
same effect. For this purpose a statistical study is made in detail of deviated variations 
of fpF2 (4foF'2) and hpF2 (dhpF 2). : 

2. Statistical Results of the Auroral F'2 Disturbances 
The ionospheric data used here are those at Troms6 (geographic lat. 69.4, geomagne- 
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tic lat. 67.1) Mar. 1956-Feb. 1957 and Jan.-Dec. 1952 and those at College (geograph. 
lat. 64.9, geomag. lat. 64.5), Mar. 1956-Feb. 1957 and Jul. 1941-Jun. 1946. As will be 
mentioned below, the variation of the f,F2 is greater and more conspicuous in the 
period of larger sunspot number. Therefore it seems to be appropriate to use the data 
during the years 1956-1957 in order to know the average daily variations of the 4f,F2 
and 4h,F2 in two stations. 
(a) Daily variations of 4f,F'2 and 4h,F2 and the seasonal change 

(1) Daily variation of 4fF2 

The daily variations of the 4{,F2 in three seasons are shown in Fig. 1 (a), (b) by 
the dotted line. It is clear from these figures that most remarkable feature of the 
variation is maximum depression of the #,F2 near noon in winter and equinox (some- 
times second maximum near 18h) and a depression throughout a day in summer. The 
hours of the negative 4f.F2 in a day become larger in summer than in equinox and 
in equinox than in winter, while the hours of the positive 4/F2 become shorter. 

We consider that these phenomena depend on the hours of the sunshine of a day 
and on the height distribution of the decay coefficient (see 4 and 5). 
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(a) (b) 
Fig. 1 (a) (b) Mean daily variations of dfF2 (dotted line, right ordinate) 
and AhpF 2 (full line, left ordinate) in three seasons. 


(ii) Daily variation of 4hp,F2 : 
For the deviation of the 4p»F2, it must be noted that the effect of the retardation 
of the radio wave by the F1 region is in some cases included in the h,F2. That is to 
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say, as the value of the f,F2 approaches that of the (,F'1 owing to geomagnetic storm 
effect, the #’F2 becomes larger and larger. This is case also for the h»F2, because 
the hpF2 is calculated by mathematical treatment for one parabolic region. In order 
to remove this retardation effect we use the same method as in the paper I. 

Process of removal of the retardation effect is applied mainly to the data in summer 
and partially those in equinox, whence the f,F2 on the quiet day is comparatively low. 
In these seasons, the retardation effect is also included in the hp,F2 on the quiet day, 
so that we must use, as a median value, the h,F2 freed from this effect, obtained by 
the process which is the same as 
that on the disturbed day. 

The daily variation of the h,- ‘fj 
F2 obtained after the above reat- 


ment is shown in Fig. 1 (a), (b) 10 
by full line. The main feature of sity is 
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season and is largest in summer 
and shortest in winter. As shown 


in the latter section, earlier maxi- 
mum of a day seems to occur by 


the superposition of the upward 
Negati ve Disturbance 


motion of the electron mass and 


the electron production at and 


after sunrise and the latter maxi- ern 


02 4 6 8 10 /2 14 16 18 2 22 % 

: : Local ti 

mum by a direction change of the Sater ae § ag ond 
Fig. 2. Mean daily variations of 4f,F2 (dotted line, 


electron mass from upwards to right ordinate) and Ah'F2 (full line, left ordinate) in 
downwards. For comparison the four stations in middle and low latitudes for two types 
daily variations of the 4fF2 and of disturbances. Data; 1949-1956 for Washington and 


’ , Kokubunji, 1938-1944 for Wathero and Huancayo. (after 
4h’ F2 in lower latitudes are shown Sato [9]) 


in Fig. 2. 
(b) Relations between 4f,f'2 and geomagnetic activity, sunspot number 

Since the magnitude of the 4fF2 is considered to depend on both geomagnetic 
activity and sunspot number (in other words, f,F'2, which is proportional to the number), 
it is better to put in statistics under a certain activity or number. Similar statistics made 
in II for Washington and other stations are insufficient, because the above consideration 
has not been taken. Fig. 3 (a), (b) and Fig. 4 (a), (b) show the mean values (4f/)F'2) of the 
Af)F2 during three hours centered on noon. In the former figure the values of the. 
AfF2, in 1942 (near the minimun sunspot number) and in 1956 (near the maximum 
one), are shown against the value of SK, (a day sum of K, indices) and in the latter 
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the Vv - > AGP, bx at S08 * F ie aay . 
alue of the 4/,F2 is shown against the Zurich provisional sunspot number with a 


different sign for different geomagnetic activities. From these figures it is clear that 
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Fig. 4 (a), (b) Relation between J4f)f2 and 
sunspot number 


(i) the negative value of 4fF2 increases with the increasing YK, in each season, 
(ii) the negative value of 4/,F2 increases with the increasing sunspot number in 
each season, and 

(iii) the negative value of 4/,F2 is largest in winter and smallest in summer. 

From (ii) and (iii) it seems that severer disturbances occur in winter than in other 
seasons and in a year of a greater sunspot number. But the above assertion is not 
true, because the percentage of decrease of the electron density seems to be about the 
same, regardless of the season and the sunspot cycle. For example, the maximum 
decrease of the f,F2 at sunspot number near 200, 120 and 40 in winter at College is 
about 5mc/s, 4mc/s and 3mc/s respectively, whereas the median values of the /.F2 
at noon on the quiet day for the corresponding numbers are 11 mc/s, 8 mc/s and 6 mc/s. 
Hence the rates of the decrease are about the same. Similar situations hold for the 
rates among the seasons. Therefore the phenomena (ii) and (iii) are ascribed to a fact 
that the fF 2 at noon in winter (especially, Nov., Feb.) is greatest in winter and smallest 
in summer and that the /,F2 increases linearly with the increasing sunspot number. 
In a year of very small sunspot number, the /,F2 is small and approaches the /F1, 
hence the decrease of the f,F2 on the disturbed day seems apparently to be small, for 
the decrease below the fl is not known. The phenomenon (i) is considered to be 
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an increase of the vertical drift velocity (see 4, 5). 
3. Theoretical Considerations 


We assume that the vertical drift of an electron on-the disturbed day is caused by 
an electric field deduced from the geomagnetic. Ds variation. When an anisotropic 
electrical conductivity is taken into account and the winds in the F and E£ regions are 
assumed to have the same directions and magnitudes, the electric current, Z, and ly, 


in the directions x (south) and y (east) are given by 


L 2 krEx “le Rugliy, 
Ty= Ryka Ry Ey; ( 1 ) 


where. oS ( Oe dz+ {ee dz CL (25 


and o,, etc. are the elements of the electrical conductivity tensor {o}= Us ad (we 
Oyx yy 

use this notation instead of {o’} in I and II). The integration in (4) is carried out 

throughout the F and E regions. E, and Ey are the electric fields in the x and y 


directions. From (1) we obtain 
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are not taken into account, the variation of the electron density is given by 


On : On 
= q(z,t)— B(z)n—v(t 
a1 q(z,t)— Blz)n Aare 


2 & 


G@) 


where ” represents the electron density, v the vertical velocity of the electron, g and 
8 the production and the effective attachment rates of the electron respectively, ¢ the 
time and z the height in unit of the scale height. ¢ is taken as Qo exp (1—z—e7*secx), 
where g is a constant and x is the solar zenith angle. The equation (7) is solved by 
the numerical method firstly attempted by Milllngton [10]. Then the equation is rep- 
laced by the following forms 


ee ESTE bw: {4 = 137 x 10% (8) 


where £=1.37x10'¢, 1/o)=1.37 10489, 1/%=Po/Qo, v=n/M, B=Bo|&\, (Bo is constant). 
If the velocity v equals to vg the solution of the equation represents the daily variation 
on a quiet day and if v equals to (vg+v,) the solution represents that on a disturbed 
day. When the former variation is subtracted from the latter the residuum equals to 


the deviation from the normal variation. 


4, Observed and Calculated Results for Individual /'2 Disturbances 


In this section we study whether the individual state of the F2 disturbance can be 
accounted for as an effect of the vertical drift of the electron. For the deduction of 
the drift velocity on a disturbed day the following points are taken into consideration. 
At first, in order to obtain the necessary Ds-variation from one geomagnetic storm the 
convenient method as shown in II is adopted. That is, the values of the horizontal 
intensity (H) and the declination (D) of the geomagnetic field for the local time nearest 
to that of the commencement of a geomagnetic storm, are taken, and from these values 
mean values of H and D on the five international quiet days corresponding to that 
hour, are respectively subtracted. The values thus obtained are the deviative parts 
at zero hour of the storm time. The process is followed for the succeeding 60 hours. 
They are denoted as 4H and 4D. Next the overlapped mean of 24 hours of them are 
calculated for the former 48 hours, assuming that the deviative values befor the com- 
mencement of the storm are zero. Then the variations can be said to correspond to 
Dst variation of an average state. When these variations are subtracted respectively 
from 4H and JAD, the parts corresponding to the Ds variation of an average state can 
be derived. aS 

Secondly we must estimate the magnitudes of the conductivity K.y and Ky. It is 
probable that the electrical conductivity of the ionosphere increases on a disturbed day 
in the auroral latitudes, owing to an increase of the ionization by the particle or other 
ionizing agencies [11], [12]. However, the height-integrated conductivity depends on the 
height interval and also an estimated height of the maximum electrical conductivity 
of the ionosphere seems to lie at or above the & region [13], [14], [15], [16]. So that 
it is unlikely that the height-integrated electrical conductivity of the ionosphere increases 
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on the disturbed day far greater than that on the quiet day, even if the electron density 
increases so much in the E or below the E region. Hence we assume that Kuy and Ky, 
are respectively twice the noon values of them on the quiet day (daily variation is the 
same as in Table 1 in II) and they are constant throughout a day. Then 
ey 1or® e.m.u., Summer ge 10" em.u., summer (9) 
7.2 x 10-8 e.m.u., equinox. y~ 15.8 x 10-8 e.m.u., equinox. 

For the calculation of the deviative part, the important point to pay an attention is the 
height distribution of the attachment coefficient 8. In I and II the distribution A in 
Fig. 5 is used and results of the calculations show that this distribution is fairly correct. 
However, the fact that as the night leng- 5 
thens (as the season approaches winter), the 


durations and magnitude of the positive values 4 
of the 4f.F2 and 4h,F2 increases remarkably 
shows that the above distribution A is not 
satisfactory, because the use of A does not 
bring the large positive value of the 4f/,F2. 
The above statistical results show that the 
rate § decreases with increasing height; for 


example, a distribution B (further discussions 


will be stated in the section 7, ii). In the pre- 


sent caulculations two distributions are used 1e ie Be 10 
6 
and the results compared. Fig. 5 Height distribution of @. Height 0 


The individual F2 disturbance (4fF2 and corresponds to the height of Chapman’s maxi- 
Ah,F2) and the corresponding calculated mate Deis 
variations (deviation of the electron density and 4h,F2) are shown in Fig. 6. The 
data in the examples are those at Tromso, 1952, as the ionospheric data in a large 
sunspot number year which we have at hand are not accompanied by the geomagnetic 
data. (we have no geomagnetic data at College). The ionospheric data on the disturbed 
day want a record much frequently, the examples are limited to some disturbances 
which have comparatively a lot of record. Now in Fig. 6 (A) represents observed results 
and (B), (C) corresponding calculated variations using the distribution of 8, respectively 
A, B and (D) the Ds variations of the 4H and 4D. From these figures it is found 
that the calculated results are well consistent with the observations. Generally the 
caulculated positive variation of the 4h,F2 is larger and continues longer for the distri- 
bution B than A and inversely the calculated negative deviation of the electron density 
is smaller. In the example (d) in Fig. 6 where the disturbance in September is shown, 
the increase of the electron density at night is seen when the distribution B is used. 
In this example the positive value of the 4f,F2 is not shown. But, as statistical results 
show a remarkable increase of the 4f,F2 in this season, it is clear that the distribution 
of @ such as B is needed. 


5. Observed and Calculated Results for Average State of F2 Disturbances 
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Since in an individual disturbance the record is few, it is difficult to compare fully 
calculated results with observation. In order to remove this defect, the variation of an 
average state and the corresponding calculated variation are compared. For this pur- 
pose, usually, a lot of calculated variations for individual state are put in statistics and 
the results compared with the corresponding statistical results of the observed data. 
But this method is much laborious and besides, the method of the derivation of the Ds 
part of the geomagnetic field, is of approximate. Therefore another process is adoped 
here. At first 4H and 4D are obtained for many, storms, they are arranged according 
to the storm time and the mean of them respectively calculated. These mean values 
are subtracted from 4H and 4D and the residuum are rearranged according to the 
local time. Then the means of them represent the Ds variations. Using these Ds 
variations and the same electrical conductivity as that in the previous section we can 
obtain a mean vertical velocity of the electron drift. Then calculated deviations during 
48 hours are obtained for the moving mass with the above velocity but with a different 
begining time. They are arranged according to the local time and the mean of them 
can be regarded as an average state. For the sake of simplicity, the representative 
beginning time is taken as 0,6, 12 and 18h of the local time. The data used for the 
calculations are those at Tromso, 1952. The calculations are made for an average state 
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AF (A | 
(A 
Otserved KME poserued 
807 DhpF2 
60 
5 if 
20- 2h [’ 
Qo 
\ 
-20- = 7 Nipeereits \A 
=~ aa Vite 
(B) 33 fof 
Calculated 


(a) Tromso, 20h, Aug. 9., 1952 (b) Tromsé, 9h, Aug. 18, 1962 
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te Iculated 


‘Geuloted 


24 hrs 
Storm time Storm time 


(c) Tromsé, 13h, Aug. 29, 1952 (d) Troms6, 20h, Sep. 25, 1952 
Fig. 6 (a)-(d) (A) represents the observed variations of 4f,F 2 (dotted line, right ordinate) and AhpF 2 
(full line, left ordinate) for the F'2 disturbance which begins at the time shown below. The time | 
is the nearest local time before the beginning of the geomagnetic storm and this is taken as zero of 
the storm time. (B) represents the calculated variation of the deviation of the maximum electro 
density* (dotted line, right ordinate) and the height (full line) using the distaibution A in Fig. 5 
(C) the calculated variation of them using the distribution B and (D) the observed variations - of 
; ae (full ie, left Gan and AD (dotted Has in y, left ee seat te a 


te hacia 
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in equinox and summer and not for that in winter. The reason is the following. If 
the ionization in the F2 region is caused by sun’s radiation as shown by Chapman, 
the ionization in high latitudes such as auroral or polar is to be zero in December. 
Nevertheless the electron density is fairly large in these latitudes. This seems to show 
that production or penetration of the electron by any other mechanisms becomes effec- 
tive in winter in these latitudes, and their mechanisms are not yet revealed. 

Average states and the corresponding calculated variations are shown in Fig. 7. 
(a), (b). In these figures (A) to (C) correspond to those in Fig. 6. Observed results 
are well analogous to those in Fig. 1 and the calculated variations are consistent with 
observed ones. It is noted that an increase of the electron density at night is shown 
in equinox when the distribution B is used, but judging from the magnitudes of the 
variations it appears that the necessary distribution of 8 is an intermediate one between 
A and B but closer to A. 


6. Physical Interpretation of Disturbance Variation 


Based on the above calculations, we state physically as for the deviations of the 
foF2 and h,F2 in the disturbance. The vertical velocity of the drift deduced from the 
geomagnetic Ds variation is upwards about from 6h to 18h and downwards from 18h 
to 6h, regardless of the season. In the absence of the electron production the upward 
motion of the electron mass with the maximum density shows a deviation in height 
which equals to the time-integration of the velocity, but if the production begins at 
sunrise the height of the maximum density rises abruptly, owing -to the production 
of an ionization layer at great height with a larger density than the maximum one of 
the electron mass which already existed before sunrise. Soon after an accumulation 
of the electron sets up below the height of the second maximum, the maximum density 
to be measured moves down to the third maximum. As an electron mass of the third 
maximum density still moves upwards, a fourth maximum density sets up below the 
height of the third, and so on. When sunset approaches an accumulation of the electron 
in the lower region becomes small, the upward drift becomes effective and the height 
of the maximum again increases remarkably. Thus the height of a maximum density 
does not so much move upwards in strong radiation of the sun (near noon), except a 
large drift velocity, and begins to move upwards, till a downward motion begins, when 
the sun’s radiation becomes weak. Therefore it is easily recognized that an interval 
between two maximums of increase in the ,F 2 is shorter in winter than in other seasons. 
In a downward motion of an electron mass with the maximum density, the height of 
the maximum falls downwards, and finally, the height remains stationary at a level 
where an abrupt increase of the electron decay begins. This is because the maximum 
electron density in the region of a large decay coefficient decreases more rapidly com- 


Ss 


pared with a decay in higher regions. : 

- On the other hand, the maximum electron density on the sunlit time (upward drift . 
in winter and equinox) becomes smaller than the median value, because an accumulation 
is difficult to take place due to the movement. The decrease of density is amplified by 
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an increase of the drift velocity. In summer when the sun is present throughout a 
day the electron density decreases compared with the median value, because of the 
penetration of the electron mass to the region of the large decay coeficient (18h-6h), 
as well as, of the difficulty of the electron accumulation. This circumstance does not 
vary when the distribution B of @ is used’ In winter the situation during the sunlit 
time is the same as that in summer, however the production is shut out on the way 
of the upward motion and the maximum electron density lying at the higher altitude 
than the median height decreases much less than the decrease of the median maximum 
density owing to the smaller decay coefficient, when the distribution B is used. So 
that an observed value on a disturbed day increases more than that on a quiet day. 
The increase continues until the density subjects to a rapid decay in the lower region 
due to a downward motion or the sun’s radiation begins to produce an ionization. 
Thus the larger the hours of the night, the larger the increase of the electron density. 


7. Discussions 


(i) It is considered that on a quiet day the height-integrated electrical conduc- 
tivity in the auroral latitudes is the same order of magnitude as that in middle latitudes, 
inferring from the geomagnetic S, variation in those latitudes and that on a disturbed 
day the electrical conductivity on the former latitudes increases whereas that on the 
latter does not vary. In the present calculations we assume the ratio of increase of 
the electrical conductivity in the auroral latitudes to be 2. This order of magnitude 
of the increase seems to be slightly small, compared with the values assumed by other 
authors. For example, the ratio is taken as 10 by Maeda [7] 10 by Fukushima and 
Oguchi [7], 3-30 by Matsushita [18] and 100 by Jacobs and Obayashi [19]. A high ratio 
is necessary to obtain enough intensity of the current to explain the observed results 
in the dynamo theory of the S, current system. But if the F2 disturbance in the 
auroral latitudes is ascribed to the effect of an electron drift as shown in this paper, 
the ratio should be in a range of 2-4 and cannot be any higher value, since otherwise 
the drift velocity on a disturbed day becomes far smaller than that on a quiet day, 
and the effect of the former drift becomes negligible. 

It may be considered.that an ionization in the aurora is fairly large and contribute 
to an increase of the electical conductivity. Indeed, the electron density shown by 
Curie et al [20] and Bullough and Kaiser [21] respectively is very large, but as shown 
by Omholt [22] the electron density in the aurora appears to be about (2-10)x 10°/cc 
and high values by Currie and others may be the ones occurred partially or temporarily. ~ 
A sporadic slant E and a blackout in the auroral latitudes are associated strongly with 
the geomagnetic bay or the bay-type variation during a storm [5] [12] and these phe- 
nomena~imply an increase of the ionization. Thus the electron density may increase 
by a few times or there about. But the level seems to lie below the level of the maximum 
electrical conductivity and the height range is probably small, hence it cannot be said 
that the height-integrated value increases with the same rate. 

Martyn states that the polarization field sets up originally in the auroral latitudes, 
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spreads over the earth and the current system due to this field has the Strength and 
the form required to explain the main features of the geomagmetic Sp, field. We do 
not then need to assume a different S, wind system from the S, and a great increase 
of the electrical conductivity necessafy to explain the intensity in the dynamo theory 
of the Sp current system. Again, a,strong electric field causes a drift of the electron 
and the F2 disurbance results which we show quantitatively in this paper. 

(ii) The increase of the (AF 2 at night on a disturbed day may be derived by taking 
the height gradient of the vertical velocity of the electron into consideration. In this 
case it is necessary that the velocity increases upwards with a great gradient. But the 
drift velocity derived from the geomagnetic Ds variation is almost independent of the 
season and hence it seems to be difficult to account for the seasonal feature of the 
increase shown by the statistics. Furthermore if the velocity gradient is present in the 
auroral latitudes it must be present in middle latithdes, also. Then it becomes difficult 
to explain the negative disturbance in middle latitudes (similar to the disturbance in 
the auroral latitudes) since the phase of the drift velocity in middle latitudes is different 
by about 180°from that in the auroral latitudes. 

The adoption of the distribution of @ decreasing upwards shown in the secttion 4 
removes the above inconsistency. This type of the distribution is analogous to the 
type suggested by Ratcliffe et a/ [23] and makes the calculated results for the distur- 
bance in middle latitudes in II better. 

We, however, do not at all abondon the presence of the height-gradient of the 
vertical velocity. On a disturbed day the gradient may be small for the large polari- 
zation field set up in the ionosphere, but on the quiet day the efiect of the gradient 
becomes remarkable. Anomalous quiet daily variation and the seasonal change of the 
electron density in the F2 region may be explained by considering the height gradient 
of the vertical velocity. Martyn [24] treats this problem under the night time condition, 
so that it remains to study the effect, including the electron production. 


8. Conclusion 


Our statistical results show that the deviation of the f.F2 in the auroral latitudes 
on the disturbed day is similar to the negative disturbance in middle latitudes, while 
the deviation of the z,F2 is different from above. The calculations of the effect of the 
vertical drift of the electron, show that the main features of the F2 disturbances in 
the auroral latitudes can be accounted for by this effect for an individual or average 
states. 

We are now proceeding with the study of the F2 disturbances ‘in the polar cap 
although the lack of the hourly values of the ionospheric and the geomagnetic obser- 
vations make it finally difficult to go on with it. 
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Abstract 


Horizontal wind systems in the ionospheric E region are deduced along the 
same line as in a previous paper [1]. This time, however, the wind systems on 
both northern and southern hemispheres in the solstice season are obtained. The 
wind systems thus obtained represent the actual state instead of the mean state 
as in Part I and II [1] [3], in the solstice season. In Part IV the effect of the 
Coriolis force is taken into account and an advance is made in the treatment in 
Part III where the wind velocity is simply assumed to be irrotational [4]. 

It is shown, as in Part III, that in summer the diurnal wind motion is 
predominant over the semi-diurnal, whereas in winter the velocity of the former 
is almost the same as that of the latter. Further the diurnal wind velocity 
is greater in summer than in winter and the semi-diurna! wind velocity is some- 
what smaller in summer than in winter. 

Some differences are found between the results of Part II] and IV. However, 
more observational material than available at present must be accumulated before 


discussing these differences, 
1. Introduction 


In Part I [3] H. Maeda deduced from the S, variation the horizontal wind systems 
in the ionospheric E region in the mean solstice. Though the daily variability and 
the heterogeneity of the ionospheric conductivity were taken into consideration in his 
calculation, he simplified the treatment by neglecting the effect of the Coriolice force 
which plays an important role in the wind motion with such a long period as that of 
24 or 12 hours. In Part II [1] the present author made an advanced treatment by 
taking the effect of the Coriolis force into account. His calculation was based on the 
same electric field in the ionosphere as that in Part I, but some differences were found 
between the results of Part I and II. 

In both Part I and II the wind systems were deduced from the data of the S; 
variation in the mean solstice, say 5 (St W), and the ionospheric conductivity in 
the equinox [3], and hence did not correspond to any actual case. Deduction of the 
wind systems on both hemispheres in the solstice season should be far significant for 
the purpose of direct comparison with the observational material. 
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In Part III and IV the wind systems on both hemispheres in the solstice season 
have been obtained. The northern and southern hemispheres correspond to the summer 
In Part III (by H. Maeda) as in Part I a simple case was 


and winter respectively. 
In this paper the effect 


treated by neglecting the effect of the rotation of the earth. 
of the rotation of the earth has been taken into consideration. The wind velocity is 


obtained by solving the dynamical equation of the atmosphere for the rotating earth 
simultaneously with the dynamo equation. 


2. Differential Equation of the Dynamo Theory 


The fundamental equations are 
E=vx H+E®, (e.m.u.) (iy 
dv a pe 
“pp t2@x¥= grad( ) (2) 


P0 
where EF and E®) are the total and the static fields respectively, v the wind velocity, 
@ the angular velocity of the earth and p’ and , the pressure variation and the static 
density of the ionosphere respectively. H is the magnetic field of the earth. As the 
static field has a potential, 
curll#=0 or E“=—grad S, tou) 


where S is a potential function. From (1) and (3) the wind velocity is given as follows; 


u=— (Ey z ) 


H; asiné 02 /, ; 
1 0S Aa) 
POMEL G re SS) 
‘ ‘ i, ( = aoog ), 


where u, v are the southward and eastward components of v respectively, and E,, Ey 4 


the southward and eastward components of £ respectively. 6 is colatitude and 2 long- 


iude. H, is the vertically upward components of the geomagnetic field which is assumed 


-to be H,=—Ccosé where ca3r. a pte the time factor as e7%, operation 
of curl on (2) gives . eh 3h Bae. | 
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Ez= 3} Een =>} {a,(0) cosm(ot+d)+b,,(0)sin m(wt +A)}, 
™m Mp 

Ey = 3) Eym= 33 (C,(0) cos m (wt +)+d,,(0) sinm (wt +A)}, 

and 


S=a3)S.=3ie¢ {A,,(0) cos m (wt +)+B,,(0) sin(w@t+a)}, 


Wm m 


then (6) is divided into the following two equations. 


aA 1 dA m? 

am ai <n ” can A An An= m 6 ’ U 
de? sin@ cos@ dé sin’0 fu) @) 

der Us : GB te ote FF =£ (4), ( 8 ) 


de? sin@ cos@ dé sin? 


where 
foe { da.» fe 1+ 2cos?0 ue 2 cos 0 (2 2 m* +2 cos*d ‘ ) 
do sin @ cos@ m do 2sin@ cosé , 
gc= -{ db,, a i+ 2.cos?6 ce 200808 ( AC» m* + 2 cos? on Jf 
dé sing cos0é m de 2sin@ cos@ 


The right-hand side of (6), or f,, and g,,, is the function of the total electric field which 
was calculated from the data of S, variations by H. Maeda [4]. 


3. Numerical method of integration 


The differential equations (7) and (8) have regular singular points at 0=0, + and 
z. In order to obtain the analytic solution over the entire range of 0, some considera- 
tions should be given to A,,, B,,, Ex,n, and Ey,,, around these singular points. 
1) around 6é=0 

We expand the solution of (7) and (8), dn, Bm, Cm, and d, as a power series in 0 about 
G=0:; 


AS = Ys boom se Ainy18 = A nyo? Sr De COLON , } 
Bn=Brm0+t Bmy8 + Bin??? + ea to ’ 

An = 4nd ate Any 9 =f Amy? Si, NAO NO ’ 

Dn = ms + Bin) + Dn ,28? + Cee Ck , (10) 
Cn = Cet CaO FC? + ges berth , 

dy, = dingo + din 510 75 Any” steer cg le tenet ae ’ 


provided that the expansion is possible. In the above formulae, A»,;, Bisjs Gris Omsir 
Cmj, and d,,,;, are coefficients of expansion with 67, Since the eastward component of 
the static field is proportional to ee and finite at 0=0, then A,,,, and By», must 
vanish. Further @,9=2,9=C2,0=d2,0=0 at @=0 according to the extrapolation of the 
geomagnetic data from the middle latitudes [3]. Substitution of (9) and (10) into (7) 
and (8) leads to 

3 Ajg=lim[30-"(@y0+ 410) + (4 d11 +5411), 

3 Byy,=lim[30-"Oy0~ Coa) + 401-54) mae a) 

lim 3 A,,0-1= —lim4 (A254 +51), 

lim 3,B95:07 = —lim A (do,1— C251), 
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In order that both sides of (11) are equal, the following conditions are given: 


Forvd;,; and £y4 


Dy,0=C1505 Q1,0= —dy,0; (12) 
for Ey,. and Ey,» ’ | 
Qo = —Aoy1, bo = Costs (13) 
and for’ Az, and 6,, 
i 1 
Aj,.= 3 fi(0) = 3 (Adi +5d;,1), 
il il 1 
By, = 3 8100) = 3 (4011-5 C141), = (14) 
Ag = 0, 
‘Rae 0, 
— together with 
Aang = myo = 0. m=1,2 


2) around $=5 


We put 0’ =5-6 and expand An, Bry Gms Ons Cm ANd dp as in 1). 
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(a) (b) 
Fig. 2 Wind systems. Northern and southern hemispheres correspond to summer and winter 
respectively. (a) Diurnal component Arrow scale: 50 m/sec for a length equal to the side of the 
small squares. (b) Semi-diurnal component Arrow scale: 25 m/sec for a length equal to the side 
of the small squares. 


Then 
a he — A’ 6 + Al ins + Ale 02 +. aris porte: eneirer : 


Bn=B' mot B mF +B nf? + Sp Ne eas ’ 


an = @' nyo eg eo ++ 556" eee bis Ti ‘ 
bn = PES BR ET PO he Si luo or ’ 
CppZC ch FC mail 10 suelo Facets ; 
Ein = mp9 AE ny F + rggQd 2 ee ee : 


(15) 


If we substitute (15) into the differential equation and equate the coefficient of 67}, 
the following condition will be given as 6>0: 


Any = —Gn50) a 
Bee hos } M= ye (16) 
3) around @=7 
The conditions as (12), (13), and (14) are readily derived by putting 6’=z—0. 


The above conditions (12) and (13) are approximately satisfied by our total electric 
field which is obtained by extrapolating the magnetic potential at middle latitudes and 


using Maeda’s conductivity [4]. In the calculation we assume that (12) and (13) are 
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Fig. 3 Distribution of pee which is proportional to the pressure variation 
wW 


Unit: 10? cm sec! 
(a) Diurnal component corresponding to Fig. 2 (a) 
(b) Semi-diurnal component corresponding to Fig. 2 (b) 
strictly satisfied. M.Hirono and T. Kitamura assumed that S, current vanishes at the 
pole [5]. Their assumption seems to be mathematically rather restrictive as known from 
(12) and (13). We have obtained a particular integral by Runge-Kutta’s method, starting 
from #0=0 and z up to ee 
The general solution of (7) and (8) consists of this particular integral and the com- 
plementary solution, viz. 
S=S,+CSp, 
where C is an integration constant, S, a particular integral and S, the complementary 
integral which is 
Sy=sinno F (7, 5 |L-+m| sin’ ) 
where F stands for the hypergeometric function. The boundary condition to be satisfied 
at the equator is that the induced field must vanish there and hence the static field 
determined by (7) and (8) is equal to the total electric field deduced from Sq Variations. 
From £, we-can determine C.* A,, and B,, for m=1 are shown in Fig. 1. 


x. : : . 
‘* From £, at omy it appears possible to determine C, However, or is necessarily qual to EF, at 
a 


ds 9 as proved in 2), because the derivative of the complementary function So vanishes at @ =F 
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4. Wind Systems 


Substraction of the static field from the total field gives the induced field, from 
which the wind velocity are readily obtained by using (4) in section 2. The wind 
systems of the diurnal and the semi-diurnal components are shown in Fig. 2. Pressure 
variatons corresponding to thesé wind systems are shown in Fig. 3. 

From these figures it is known that in summer the diurnal wind velocity is greater 
than the semi-diurnal, and in winter both diurnal and semi-diurnal wind velocities have 
almost the same magnitude. In Part I and II we obtained the wind systems in the 
mean solstice. The diurnal wind velocity was shown to be greater than the semi- 
diurnal. Therefore, actually, the mean solstice corresponds rather to summer. 

If the wind velocities of both solstice seasons are compared, it is found that the 
diurnal wind velocity is greater in summer than in winter, whereas the semi-diurnal 
wind velocity is a little greater in winter than in summer. The diurnal velocity is 30~ 
40 m/s in summer and 30~15 m/s in winter in middle latitudes. The semi-diurnal 
velocity is 20~30 m/s in winter and 10~15 m/s in summer in middle latitudes. As in 
Part II, the cyclonic character is noticeable in the diurnal component. The semi-diurnal 
pressure variation reaches the maximum value earlier in the morning in summer than 


in winter. 


5 Discussions 


In recent years the observation of the wind motion in the ionosphere has been made 
by using various techniques at some locations [6] [7] [8] [9]. However, consistent results 
have not been attained as yet, presumably because of difficulty in the analysis of the 
data. 

Our calculation shows that in summer the diurnal wind velocity is greater than the 
semi-diurnal ; the former being about 2 times greater than the latter in middle latitudes. 
Further it is shown that the diurnal wind velocity is greater in summer than in winter. 
These results seems to roughly agree with those obtained from the wind observation 
by L. Harang and K. Pedersen [9]. However, our results are discrepant with what we 
expect from observation of the daily pressure variation on the ground as well as with 
an observational result obtained by J. S. Greenhow and E. L. Neufeld in England [8]. 
Their results showed that no regular diurnal velocity had been detected. 

The ionospheric conductivity used in our calculation is based on the assumption 
that the electron density at night is 1/12 of that at noon. If this value is unavailable, 
our calculation should be revised. Nocturnal ionization has recently been studied by 
A. P. Mitra on a semi-empirical basis [10]. His estimate seems to agree with that value. 

There is a considerably rapid shift of the phase with height in the atmospheric 
oscillation at the bottom of E region. Hence, as suggested previously [1], our result 
obtained by taking mean with height may not represent the actual state. However, a 
simple calculation shows that any conspicious change will not come in our present result 
when the shift of the phase is taken into consideration [11]. 
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It has been assumed that the wind motion in the ionosphere is altogether periodic. 
If non-periodic motion exists, our result will suffer some alterations. 

At present the results of various wind observations seem somewhat imcompatible 
with each other. It is desirable that wind observation be conducted at eny places. 
Further theoretical study is very important. 
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. Addendum 


Before receiving the first galley proof the results of a systematic investigation of 
winds in the upper atmosphere had been reported (W. G. Elford, Mixed Commission 
on the Ionosphere, New York, August, 1957). According to that result the periodic com- 
ponents contain the diurnal as well’ as the semi-diurnal terms. Furthermore, the order 
of magnitude, the seasonal variation, and the sense of rotation, of the wind vectors 
are roughly agreeable with the result of our calculation. 
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LETTERS TO THE EDITORS 


Distribution of Ions around ‘charged Fine Wires 


According to the results of the writer’s measurement* of the dissipation of electricity 
from charged fine wires, the dissipation coefficients are given by 


a=a(1+2), (1) 


where ‘a’ is the dissipation coefficient, a, 6 are constants, and 7 the radius of wire. 
On the basis of this relation some considerations on the distribution of ions around the 
charged fine wires were made. 

When the wire is charged positive the negative ions would be collected around it 
due to the intense electric field, and the positive ions repelled away. The charges 
entered into the wire of unit length in time interval dt will be given by 


dq’ =2nrpkE dt, P hey 


where dg'=charges flowing into the wire, a 
o=density of charge at the surface of the wire, 
E=intensity of the electric field, 
k=mobility of ions. 
Since, 
rE=2q=2CV, | oe ee 


where q is the charge of the wire per unit eens at any time, C is the capacity of it 
_ with respect to the earth and V is the potential, from (2) and (3) : 
eee ee estab dal =Ankpqdt. “8 a 
On oa ide the charge on the wire dissipates dg in df, then 
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wire is given by : 
1 


C=—__————. for h)r, 8 
2log, 2h/r po ey 


where h/ is the distance between the wire and the earth. Within the limits of values 
of r (0.044~0.0063cm), which were used in writer’s measurements, C is briefly represented 
by the next approximate formula; 


C=a+6r, ’ (9) 
where a and £ are constants, whose values are 0.073 and 0.68 respectively. From (7) 
and (9) we obtain 

p=no(1+5), | | (10) 


- neglecting small quantities. Here po) and » are as follows; 


0o= bo +4 p/aBV, = (11) 
eet ae (12) 
(0. +4y'aBV 


These constants should be determined in connection with the measurement. Combinning 
(6) and (10), writing 4zkp)>=a, we get | 


_— 
—- 


| Ss San obi a= Ankpo(1+2) 2 ES 


— 4 r 4 ei ay 


S| <a tua erate nn acti Ce) 


Taking »=6, equation (13) is identical with (1). The term 4/a8V in (11) is small com: 
- pared with o/, therefore (11) and (12) becomes 


= as 
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“In (14) it is cles 
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or 
n=n(1 +2), (i=). (17) 


where ,=a/4zke and e is the ele- 
mentary charge respectvely. By 
(17) computations of ion number 
around charged wires are pos- 
sible. In Fig. 1, it is shown the 
ion number for k=1 and 0.5 (cm/ 
sec/volt/cm), taking as the mean 
value a,=2.00x10-* (1/sec) (= 
2.00x10-*/42k) and b6=0.0074 
from the measurement. For large 
Yr, b/r->0, therefore n= Oo/e) in 
_ (17) may be regarded to show the 
number of ions around the con- 
ductor of large r. On account 
of this the 7) have been taken 
as the ion number i in free atmos- 
phere. And it is noteworthy that the values in Fig. 1 are the one “Obtained in the 
laboratory without air circulation. 7 
The writer is indebted to Professor ae Tamura of Kyoto University for his i in- a 
valuable advice and EOIN {- | 
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Wind Systems for the Geomagnetic Sa Field 


4 


Horizontal wind-systems in the ionospheric E-region required to produce the solar 
daily magnetic variation on disturbed days (say, S») are deduced along the same line as 
in the case of S, [1, 2]. The data 


used are the horizontal intensity of 120 a oX-COMA 
; °Y come 


the geomagnetic S,-variations* at $d, 
62 stations over the world for the 2 
mean solstice during the Inter- 
national Polar Year, 1932-33. The 
harmonic coefficients of the North 
(4X) and East (4Y) components are 
plotted in Fig. 1. The dotted (for 
4X) and full (for 4Y) lines show 
a smooth distribution of the coef- 


ficients, where the curves between 
60°N and 60°S are derived from 
the potential of the S, [3] and Sp 
[4] fields* for the mean solstice (see 
Table 1), and the curves in latitudes 
higher than 60° are deduced direct- 
ly from the distribution of the 
coefficients. 

The distribution of electrical 
conductivity on disturbed days is 
assumed as follows: 


. Dee baer (¢) 
yy t=} Ly (Y) ¥(¢,t), 
Dipy xy0 Y 


where ¢ is the latitude, and ¢ the 
local time. The values of >},’s are 
given in Table III (Model B) of a 


* In this paper, the disturbance-daily 
variation Sp is defined as the average N GEOMAGNETIC LATITUDE 


- daily vouaeo ne ave Pe eee : Fig. 1 Latitude-distributions of the harmonic coeffici- 
disturbed minus quiet days (i. €., Sa ents for the geomagnetic North (4X) and East (AY) 
Sq), so that Sg=Sg+Sp, and Sp components of the geomagnetic S,-variations for the 
seems an average of the diurnally mean solstice. The dotted (for 4X) and full (for AY) 
varying part D, of magnetic storms. lines show a smooth distribution of the coefficients. 
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Table 1. Spherical harmonic coefficients for the geomagnetic Sg- and Sp-fields 
analyzed in the form 
V=nH DD (am cos mt+bm sinmt) Pr (cos A), 
nm m 


using the Second Polar Year data. (Unit=7) 


Sq Sp 
m | n ax ba n ax bn : 
oss = = ; 
2 8.40 —2.35 2 -0.71 =6.55 . 
1 4 His 0.40 4 —0.23 —1.97 
6 0.21 0.12 6 0.53 ~0.46 
3 ete 1.69 3 2.32 —3.46 
oN hog 0.14 =0.15 5 ~0.94 0.38 
7 —0.02 0.22 7 0.06 ~0.69 
4 1.42 ~0.86 4 ~0.38 0.57 
3 6 0.30 —0.02 6 i.ites 0.53 


previous paper [5], and the latitu- 

Ce ae Rr _ de distributions of Y(y,t) at inter- 

j I/\ bo vals of two hours are illustrated in 

: Sea Ree eS cs * on Fig. 2 which is reasonably deduced 
ee Ae : 


= . Ve from f,E and f,E, of the ionos- 
9-12" (noon) phere at sunspot-minimum. 

The total electric field for S,; 
is calculated by using the geo- 
magnetic intensity and the con- 
ductivity prepared above, and the © 
dynamo equation is solved, con- | 
sequently S; wind-systems are — 
obtained. Since the method of — 
the SprindslZk ie detailed des 
criptions are not repeated here. 
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Fig. 2 Latitude-distributions of ‘es 
of uo bburs. 
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Fig. 3 Wind-systems for the Sj-field, viewed from above the North pole: 
(A) the diurnal and (B) the semi-diurnal component. 
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Fig. 4 Wind-systems for the Sg-field, viewed from above the North pole: 
(A) the diurnal and (B) the semi-diurnal component. 
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Meeting of the Society of Terrestrial Magnetism and Electricity : 4 

The 22nd General Meeting was held at the Kakioka Magnetic Observatory oil 

October 9-11, 1957. 

Number of the Reports read at the Meeting: | 
Geomagnetism, 16; Ionosphere, 6; Night Airglow, 4; Cosmic Rays, 11; f 


Rock Magnetism, 5; Earth Current, 2; Atmospherics, 3; 


Atmospheric Electricity and Radio Meteorology, 7 
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